First order shape transition and critical point nuclei in Sm isotopes from relativistic 
mean field approach 
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The critical point nuclei in Sm isotopes, which marks the first order phase transition between 
spherical U(5) and axially deformed shapes SU(3), have been investigated in the microscopic 
quadrupole constrained relativistic mean field (RMF) model plus BCS method with all the most 
used interactions, i.e., NL1, NL3, NLSH and TM1. The calculated potential energy surfaces show 
a clear shape transition for the even-even Sm isotopes with N = 82 ~ 96 and the critical point 
nuclei are found to be 48Sm, °°Sm and 1°?Sm. Similar conclusions can also be drawn from the 
microscopic neutron and proton single particle spectra. 


PACS numbers: 21.10-k, 21.60.Jz, 21.60.Fw, 21.10.Re 


The equilibrium shape for finite many-body systems 
such as nucleus, atom, molecule, etc. have been a hot 
topic for the past several decades. These shapes are nor- 
mally not rigid and change as a function of their con- 
stituents. Sometimes these changes can be quite abrupt 
and exhibit phase transitional and critical point behavior 
similar to that found in a wide variety of many-body sys- 
tems, though the concepts of phase transitions are only 
approximate for finite system. 


Properties of the system in the transition region and, 
in particular, at the critical point can be found by solving 
numerically the eigenvalue problem for the Hamiltonian. 
However, the transition regions are difficult to interpret 
as they exhibit a complicated interplay of competing de- 
grees of freedom. Yet such system are in many respects 
the most important, as their structure defines the nature 
of the transition region itself. 


Recently, the concept of such critical point solutions 
has been introduced in the framework of the algebraic 
models, in which different shapes (phases) correspond 
to dynamic symmetries of some algebraic structure G 
and the phase transition corresponds to the breaking 
of the dynamical symmetries. In the Interacting Boson 
Model [1] G = U(6) and there are three dynamical sym- 
metries characterized by the first algebra in the chain 
originating from U(6), that is U(5), SU(3), and SO(6), 
with spherical, axially deformed, and y-unstable shapes. 
Experimental examples of all three types of symmetries 
can be found in nuclei, e.g., the U(5) nuclei can occur 
near closed shells, SU(3) nuclei can occur in the mid- 
dle of shells, and the O(6) limit tends to occur in nu- 
clei with particles-holes configuration. The phase transi- 
tion between spherical U(5) and axially deformed shapes 
SU(3) is first order [2, 3]. The phase/shape transition 
region is characterized by a pronounced 8 softness [4]. 
This characteristic led to the analytic solution for critical 


point nuclei, called X(5) for axially symmetric case, and 
is based on analytic solutions of the Schrödinger equa- 
tion corresponding to a geometric (Bohr) Hamiltonian 
with a square-well potential [5]. It has been shown that 
152Sm [6] and other N = 90 isotones [7] are the first em- 
pirical manifestation of the predicted nuclei at the critical 
point of vibrator (spherical shape) to axial rotor (axially 
deformed shape) phase transition. 


The critical-point symmetries provide a classification 
of states and analytic expressions for observables in re- 
gions where the structure changes most rapidly. The 
ideas can be visualized in terms of two coexisting poten- 
tials, spherical and deformed, whose energy separation 
varies with nucleon number. Although phase transitions 
and critical point have been discussed in atomic nuclei 
for many years, there is no microscopical investigations 
yet. As an important complementary way, a microscopic 
study could be very helpful for us to understand where 
and how nuclei exhibiting structures close to these sym- 
metries occur. It is the main purpose of this Letter to 
close this gap in the study of symmetry. The relativistic 
mean field (RMF) theory [8] has received wide attention 
because of its successful description of many nuclear phe- 
nomena during the past years [9]. In particular, the ex- 
otic halo phenomena can be understood self-consistently 
in this microscopical model after the proper treatment of 
the continuum effect [10, 11, 12]. Here in this letter we 
will discuss the realization of the critical point concept in 
atomic nuclei in the microscopic quadrupole constrained 
RMF theory with pairing treated by the BCS method, 
which will show that 148:150,152Sm are excellent empiri- 
cal manifestations of this critical point structure. 


The basic Ansatz of the RMF theory is a Lagrangian 
density where nucleons are described as Dirac particles 
which interact via the exchange of various mesons in- 
cluding the isoscalar-scalar 0 meson, the isoscalar-vector 


TABLE I: The total binding energy and the quadrupole de- 
formation 62 of the ground-states of ‘“4~'°8Sm calculated by 
the constrained RMF theory with the parameter sets NL1, 
NL3, NLSH and TM1. 


E[MeV]| EXPO NL1 NL3  NLSH TM1 
Sm | 1195.74 1201.28 119834 1200.20 1201.63 
146Sm | 1210.91 1213.83 1211.79 1214.00 1215.89 
1485m | 1225.40 1225.95 1225.48 1227.95 1230.12 
150Sm | 1239.25 1239.59 1239.29 1241.09 1243.54 
1525m | 1253.11 1252.78 1253.54 1255.39 1256.12 
1545m | 1266.94 1265.50 1267.54 1269.59 1268.99 
1565m | 1279.99 1276.98 1280.05 1282.26 1281.61 
1585m | 1291.98 1287.57 1291.73 1293.99 1293.07 


NL3 NLSH TM1 

0.00 —0.01 
46Sm 0.00 0.05 0.06 0.08 0.06 
148 Sm 0.16 0.15 0.14 0.14 0.13 
1505m 0.21 0.25 0.23 0.21 0.16 
1525m 0.24 0.34 0.30 0.28 0.19 
1549m 0.27 0.35 0.33 0.32 0.31 
156 Sm 0.28 0.36 0.34 0.33 0.33 


158 Sm 0.28 0.36 0.35 0.34 0.34 


w meson and the isovector-vector p meson. The effective 
Lagrangian density considered is written in the form: 


L = ti (iP - M) Wi + 300o — U (o) — gofio yi 
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= ZEE = bi Adi, (1) 
where y = 17° and w is the Dirac spinor. Other sym- 
bols have their usual meanings. 

The Dirac equation for the nucleons and the Klein- 
Gordon type equations for the mesons and the photon are 
given by the variational principle and can be solved by ex- 
panding the wavefunctions in terms of the eigenfunctions 
of a deformed axially symmetric harmonic-oscillator po- 
tential [13] or a Woods-Saxon potential[14]. The details 
can be also found in Ref. [9] and references therein. 

The potential energy curve can be calculated micro- 
scopically by the constrained RMF theory. The binding 
energy at certain deformation value is obtained by con- 
straining the quadrupole moment (Q2) to a given value 
jig in the expectation value of the Hamiltonian [15], 


(H) = (H) +50. ((@2)- ma), ) 


where C, is the constraint multiplier. 

For the nuclei studied in this paper, the full N = 20 de- 
formed harmonic oscillator shells are taken into account 
and the convergence of the numerical calculation on the 
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FIG. 1: The potential energy curves for "447158 Sm obtained 
by the constrained RMF theory with the NL1 parameter set. 
The ground state binding energy is taken as a reference. 
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FIG. 2: The same as Fig. 1, but with NL3 


binding energy and the deformation is very good. The 
converged deformations corresponding to different u2 are 
not sensitive to the deformation parameter ĝo of the har- 
monic oscillator basis in a reasonable range due to the 
large basis. The different choices of $9 lead to different 
iteration numbers of the self-consistent calculation and 
different computational time. But physical quantities 
such as the binding energy and the deformation change 
very little. Thus the deformation parameter Bo of the 
harmonic oscillator basis is chosen near the expected de- 
formation to obtain high accuracy and low computation 
time cost. By varying u2, the binding energy at different 
deformation can be obtained. The pairing is considered 
by the constant gap approximation (BCS) in which the 
pairing gap is taken as —= for even number nucleons. 


VA 


The binding energies and the quadropole deforma- 
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FIG. 3: The same as fig. 1, but with NLSH 


tion for the ground states are listed in Table I in con- 
strained RMF theory with NL1[16], NLSH/[17], TM1[18] 
and NL3[19] parameter sets. For the binding energies, 
the data are well reproduced within 0.5 %. Particularly 
for NL3, excellent agreement (within 1 MeV) is observed 
for the binding energy in 46-158Sm. Even for neutron 
magic nuclei 144Sm, the difference between the RMF cal- 
culation and the data is less than 3 MeV, i.e., relatively 
less than 0.3 %. The spherical shapes in '44146Sm and 
the weakly deformed !48Sm are well reproduced. The de- 
formations in °-'58Sm are a little over-estimated by the 
theoretical calculations. In general, it can be concluded 
that the data is well reproduced by the constrained RMF 
theory. 

Figures 1, 2, 3, and 4 show the potential energy 
curves for 44-158S5m in constrained RMF theory with 
NL1, NL3, NLSH and TM1 parameters, in which the en- 
ergy for the ground state is taken as a reference. Similar 
patterns are found for all the parameter sets. The ground 
state of '44Sm is found to be spherical and has about 
12 MeV stiff barrier against deformation. Although the 
ground state of '4°Sm is still spherical, its barrier be- 
comes around 8 MeV against deformation. With the in- 
crease of the neutron number, the ground state gradu- 
ally moves toward the deformed side and the potential 
energy curve becomes more soft. Finally well deformed 
154—158Sm are observed. 

To examine how the shape of the Sm isotopes changes 
with the neutron number, the differences of the bind- 
ing energy between the ground state and the state with 
spherical shape for Sm isotopes are presented in Ta- 
ble II in the constrained RMF theory with the NL1, 
NL3, NLSH and TM1 parameter sets. These differences 
can tell us how soft the nucleus is against deformation 
and may give us a hint on the phase transition of the 
nuclear shape. From !44—1°8Sm, the energy differences 
between the ground state and the state with spherical 


TABLE II: The difference of the binding energy between the 
spherical state and the ground state in unit of MeV calculated 
by the constrained RMF theory with the NL1, NL3, NLSH 


and TM1 parameter sets for 1447158 Sm. 
NL1 NL3 NLSH TM1 
Sm 0.00 0.00 0.00 0.00 
466m 0.02 0.11 0.33 0.07 
148 Sm 2.26 2.67 2.72 1.98 
1505m 6.11 5.19 4.13 3.08 
1525m 9.51 8.02 6.59 3.38 
mgm 12.33 10.73 9.20 4.37 
156 Sm 14.10 12.61 11.64 6.53 
1585m 14.80 14.36 14.51 8.68 
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FIG. 4: The same as Fig. 1, but with TM1 


shape change from 0 to 15 MeV. There are two jumps in 
the energy differences. The first jump appears at 1485m 
and the second at 154Sm, which suggests the shape tran- 
sition from spherical to critical point nuclei, and finally 
to axially deformed nuclei. The potential energy curves 
for 1485m, 1°°Sm, and 15?Sm are relatively flat, i.e., they 
are (2-soft nuclei in the transition between spherical and 
axially symmetric deformed nuclei. 

One of the merits of microscopic nuclear models such 
as RMF theory is that it can provide detailed informa- 
tion on single particle levels, shell structure etc., which 
are very important for us to discuss nuclear structure 
and examine the deformation driving effect. In Fig. 5 
the single neutron levels for 44-'8Sm lying between 0 
and —15 MeV are shown. The Fermi levels are presented 
by dashed line. Fig. 6 is similar to Fig. 5 but for pro- 
tons. These two figures present results calculated with 
the parameter set TM1. The other three parameter sets 
give similar single particle structure thus not presented 
here. 

From Figs. 5 and 6 one finds a consistency between the 
shell structure evolution and the shape evolution in Sm 
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FIG. 5: The single neutron levels for ‘“4~!°8Sm obtained by 
the constrained RMF theory with the TM1 parameter set. 
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FIG. 6: The single proton levels for ‘*4~1°8Sm obtained by 
the constrained RMF theory with the TM1 parameter set. 


isotopes with increasing neutron number. Namely, for 
449m, spherical symmetry is better restored. The de- 
formation develops with N increasing. The deformation 
in 1465m is still small and the energy gap with N = 82 
can be clearly seen from the single particle spectra. The 
nuclei 148-150Sm belong to a transition area, in which 
the N = 82 gaps still exist but much smaller than that 
in 144146Sm, Starting from 4Sm, the gap with N = 82 
disappears. Meanwhile a deformed gap for Z = 62 de- 
velops and also that for N = 94. Correspondingly, we 
observe the well deformed !°4~!°8Sm. 

In summary, the shape transition between spherical 
and axially deformed nuclei in '44~1°8Sm are investigated 
by the microscopic quadropole constrained relativistic 
mean field theory with all the most used interactions, 


i.e., NL1, NL3, NLSH and TM1. The RMF calculation 
reproduces very well the data of the binding energy and 
deformation for the ground states. By examining the 
potential energy curves and the single particle levels ob- 
tained by this microscopic approach, 1485m ,!°°Sm, and 
152Sm are found to be soft against 3 deformation and the 
critical point candidate nuclei, which marks the first or- 
der phase transition between spherical U(5) and axially 
deformed shapes SU(3). 
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